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Scanning tunneling microscopy can provide a probe for the detailed study of quasiparticle states 
in high- Tc superconductors. We propose that it can also be used to acquire specific information 
about impurity-induced quasiparticle states and the superconducting order-parameter structure. In 
particular, the local density of states is found to be sensitive to impurity-induced resonances and to 
the symmetry of the order parameter. 



There has been much recent interest in the role of im- 
perfections in unconventional superconductors, because 
they may provide a fruitful direction for understand- 
ing the pairing mechanism in high-T^ superconductors. 
Nonmagnetic impurities are known in conventional su- 
perconductors to perturb the superconducting proper- 
ties only weakly [1]. However, in superconductors with 
a nontrivial order-parameter phase, the same impurities 
are shown to be pair-breakers [2] producing a finite life- 
time for quasiparticles near the nodes and a finite density 
of states at low energies. Impurities, therefore, provide 
practical and stringent tests for theoretical paradigms de- 
scribing high-Tc superconducting materials, and the lack 
of consensus regarding the proper identification of the 
pairing state provides a clear impetus for studying impu- 
rity effects further in more general situations. 

Questions concerning spatially inhomogeneous quasi- 
particle states generated by impurities were first explored 
in the 60 's [3]. It was shown that a magnetic impurity 
in an s-wave superconductor, interacting with the spin 
density of conduction electrons, produces a bound state 
inside the energy gap. In general, the overlap with the 
particle-hole continuum only allows virtual states to be 
formed with finite lifetimes. Similarly, it has been re- 
cently emphasized that, in a rf-wave superconductor, a 
single scalar impurity generates a low-energy resonance 
which, in the unitary limit, becomes a marginally bound 
state of a highly anisotropic character [4]. 

In this Note, we study the effect of strongly-scattering 
impurities on quasiparticle states at low energies in a 
(i-wave superconductor and discuss how these modi- 
fied states could be observed in scanning tunneling mi- 
croscopy (STM) which gives access to a much richer set 
of diagnostic tools by probing the local density of states 
[5]. To this end, we present a detailed description of 
the density of states as a function of space and energy 
near an impurity. Our main results are twofold: first, 
we emphasize that a scalar impurity breaks locally the 
particle-hole symmetry leading to experimentally observ- 
able consequences and, second, we show that the reso- 



nance behavior arising from the virtual impurity-bound 
state strongly enhances the local density of states at the 
resonance energy in the vicinity of an impurity. This en- 
hancement is sensitive to the order-parameter structure, 
and its anisotropic spatial dependence refiects the orbital 
symmetry of the pairing state. Because the local density 
of states can be determined experimentally by measuring 
the local I{V) characteristics, we suggest that STM stud- 
ies can also reveal the nontrivial impurity-induced states 
in (i-wave and other unconventional superconductors. 

The local density of states of a rf-wave superconduc- 
tor near an impurity has been studied earlier both in the 
quasiclassical theory [6] and in the Born approximation 
where the impurity scattering acts as a weak perturba- 
tion [7]. It was found that impurity scattering modifies 
the local density of states in a manner which may serve as 
a finger print of rf-wave pairing. While weakly-scattering 
impurities do not produce a strong dependence on the 
phase of the pairing state, the resonance behavior of the 
virtual impurity states in the energy domain allows, for 
example, the phase-sensitive contribution of the pairing 
state to be isolated and explored. In addition to the 
Josephson effect, this provides a potentially attractive, 
new method for probing experimentally both the ampli- 
tude and the phase of the order parameter of high-T^ 
superconductors [8]. 

Specifically, we consider a repulsive (5-potential impu- 
rity in a superconductor with a real gap function A(k): 
U{v) = Ub{v), U > 0; our results can easily be general- 
ized for an attractive impurity. Scattering of quasiparti- 
cles from the impurity is described by a T-matrix, T{uj), 
which is independent of wave vector for s-wave scattering. 
The Green's function in the presence of the impurity is 
G(r, r'; uj) = G^°Hr - r', lo) + G^°Hr, Lo)f(Lo)G^°H-r',Lo), 
where both &'^\r,Lj) and T(lj) are matrices in Nambu 
particle-hole-spinor space, with w-l-iO"*" implied hereafter. 
The Green's function of a clean system is G('^)(k,w) = 
[wfo — e(k)f3 — A(k)fi]~"'", where e(k) is the quasiparticle 
energy, f„ (a = 1,2,3) are the Pauli matrices, and tq is 
the unit matrix in Nambu space. 
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For a system with rf-wave pairing and a particle-hole 
symmetric normal-state energy spectrum, the T-matrix 
has only two nonzero components [9]: T = TqTo + Tsfs. 
They are To(lo) = Go(lo)/[P - Goiojf] and Taioj) = 
c/[P — Go(u))^], where the dimensionless parameter c 
is inversely proportional to the strength of the impu- 
rity potential U: c = {i:NpU)~^. The momentum- 
averaged Green's functions, Ga(w) = G^a\v = 0,w), 
are conveniently expressed in terms of their real-space 
analogues GL°^(r,w) = /c^'^^ Tr f„G(°)(k, w)e''' 

(a = 0, . . . , 3); Np is the density of states per unit vol- 
ume at the Fermi energy in the normal state. The virtual 
or bound states in the single-impurity problem are given 
by the poles of the T-matrix: Go(Sl) = ±c. Choosing the 
gap function of dj;2_y2 symmetry at the Fermi surface so 
that A(i^) = Ao cos 2ip, we obtain 

Go(c.) = -(2e/7r) [sgn(c.)A'( ^1 - e^) + iKie)] , (1) 



for e = |w|/Ao < 1. Here, K is the complete elliptic 
integral of the first kind. In principle, the solution for 
the poles is complex, indicating the resonant nature of 
the virtual state. For small c, this is easily seen from 
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which solve SI = =fS1o — *T to logarithmic accuracy. How- 
ever, as c ^ 0, the resonance becomes arbitrarily sharp. 
For Born scattering, c ^ 1, no resonance structure is 
generated. 

STM measurements can probe the local density of 
states by measuring the differential conductance dl/dV 
as a function of bias voltage V . We focus on the den- 
sity of states at zero temperature defined as N(r,ui) = 
— Im Gii(r, r; w). (Hereafter the density of states of a 
superconductor is normalized relative to Np.) It is con- 
venient to divide it into two parts: N(r,ui) = No(uj) + 
6N(r,ui), where No(uj) = — ImGo(w) is the uniform den- 
sity of states of a clean superconductor and 6N(r,ui) 
is the impurity-induced contribution. For \ui\ ^ Aq, 
No(uj) = |w|/Ao. Using the above results, the impurity- 
induced contribution can be written as: 



6N (r,uj) 



-Im[Gf(r,c.)T+(c.)Gf (-r,c.) 

+ G[°\r,u,)T.{u,)G[°\-r,u,) 



(3) 



where we have defined Gf^ = G^ 
To ± T3. In the limit c ^ co (i.e., U ^0), this equation 
reduces to the Born form [7], where To = 0(c~^) and 
T3 = 0(c~^). However, in the opposite limit, the reso- 
nance behavior becomes an important new ingredient in 
the problem, as can be seen by noting that 
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Thus, in the neighborhood of w = ±Slo, a rapid modula- 
tion of the differential conductance is expected. Depend- 
ing on the relative magnitudes of the resonance energy 
Slo (determined by c) and the probe energy uj (i.e., the 
bias voltage), three regimes of interest appear. Here, for 
brevity, we consider only the strong-scattering limit, for 
which Slo ^ Ao. (i) For small energies \ui\ ^ Slo, the 
problem reduces to that of a weakly scattering impurity, 
described by the Born approximation [7]. In this case, 
T+ ~ — T_ . (ii) In the opposite limit, where |w| ^ Slo, we 
have T+ ::±T-. (in) In the cross-over regime, \uj\ ~ Slo, 
they are rapidly varying functions of w. In fact, T± has 
a resonance at w = =fS1o- This resonance is important 
already at a moderately small value of c. As c decreases, 
the amplitude of the resonance peak in ImT± increases 
rapidly and, at the same time, the resonance becomes 
sharper but this happens only at a much slower rate as 
is evident from Eq. (2). 

Recall the three characteristic scales of length: the 
impurity-induced length Iq = Tivp/^o, the coherence 
length ^0 = ^vp/Ao, and the Fermi length £p = 2'iT/kp, 
where vp is the Fermi velocity and kp is the Fermi wave 
vector. The applicability of the continuum theory and 
for having a pronounced impurity-induced resonance re- 
quire that £p 't^ Co 't^ io- The nature of the local density 
of states depends on the ratio of the distance r from the 
impurity to these scales. Below, we focus on N(r,ui) in 
the vicinity of the impurity because, for r > £0, the inten- 
sity of the STM signal is rapidly saturating function of 
r, and the tunneling current will be essentially spatially 
uniform. 

(i) A d-wave supereonduetor. First, consider the differ- 
ential conductance at the repulsive impurity site, r = 
(i.e., r ^ ip). Because G(°)(r = 0,w) = Go(w)fo, the 
density of states has a simple form: N(v = 0,w) = 
— Im [cGo(w)/(c — Go(w))]. For finite c, the van Hove 
singularity in the original density of states is replaced by 
zero at w = ±Ao. Overall, for a weak impurity, c ^ 1, 
the density of states is approximately the same as for a 
clean superconductor: N(v = 0,w) = #0(1^) -|- 0(c~^). 
The leading-order impurity effect is the spectral-weight 
transfer within the energy gap region (\uj\ < Ao) from 
positive to negative energies. In contrast, for a strong 
impurity, c ^ 1, the situation becomes much more in- 
teresting. For \ui\ ^ Slo, the density of states is neg- 
ligible: N(r = 0,ui) = 0(c2), whereas, for \ui\ < fio, 
we recover the unperturbed density of states. In the 
crossover regime and for uj < 0, N(r = 0,w) has a 
clear peak of magnitude approximately proportional to 
clogc""'^, whereas, for w > 0, it shows only a very shal- 
low feature; see Fig. 1. As c ^ 0, the weight of both of 
these features vanishes continuously, and their position 
scales with Slo to zero energy. The resonance peak at 
UJ < persists up to c ~ 1 where it appears as a broad 
structure. It is natural to interpret this peak as a refiec- 
tion of the impurity-induced quasiparticle states, which 
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FIG. 1. The local density of states N{r = 0, to) of a rf-wave 
superconductor for (a) a; > and (b) a; < at the impurity 
site as a function of energy to for c = 0.02, 0.06, and 0.10 (from 
bottom to top). Here, Qo = 0.03Ao is the impurity resonance 
energy for c = 0.10. The density of states No{io) = |a;|/Ao of 
a clean rf-wave superconductor is shown by the dashed line. 
Note the different density scales. 

are hole(electron)-like for a repulsive (attractive) poten- 
tial. Consequently, an STM experiment will see only a 
small differential conductance at low energies (|w| ~ Slo), 
and approximately zero conductance at higher energies. 
The depletion of the quasiparticle states is a result of a 
strong impurity potential which allows quasiparticles to 
be excited essentially only at very high energies \ui \ ^ Aq 
[10]. Because particles predominantly of hole (c > 0) or 
electron (c < 0) character can reside at the impurity site, 
a strongly-scattering impurity leads simultaneously to a 
nearly complete local pair-breaking effect. In the uni- 
tary limit (c 0) and at finite energies (|w| ^ the 
density of states vanishes for both signs of the impurity 
potential in both the electron and hole channels. 

Next, consider the density of states in the neighbor- 
hood of the impurity, Ip ^ f < £o, the focus being at 
distances r ~ and energies |w| ~ SIq ^ Aq. Thus, in 
computing the Green's functions, we may neglect contri- 
butions proportional to uj/Aq, because they become im- 
portant only at large distances r > £0 or at higher ener- 
gies [11]. Along the gap extrema, their dominant part can 

be estimated as G'[°|(r,w) ~ (£f /r)i/2e-''/^° 0(£f/j'). 
Clearly, along the diagonal directions where the energy 
gap has a line of nodes, G^'f ^ vanishes. In contrast, 
G!|^-'(r,w) ~ Co/f. This result, valid for distances r ^ 
= hvp/lcj], is the leading contribution in &'^\r,Lj), 
for |w| ~ Slo and ^0 < f < to- This type of behavior, but 
extending to arbitrary long distances and refiecting the 
impurity-induced marginally bound state [4], can be ob- 
served only in the unitary limit where we first take c ^ 
and then w ^ 0, because ImT±(0) = 0, for c 7^ 0. 

A number of important and novel results follow from 
the qualitative behavior of T± and First, it is in- 




FIG. 2. The local density of states N{r,Lo) of a rf-wave 
superconductor for (a) a; = Qo and (b) to = — Qo as a 
function of position r = {x,y) around a repulsive impurity 
(c = 0.1) placed at r = 0. The impurity resonance energy 
is Qo = O.O3A0 and S,o/iF = 6. The maximum values of 
7V(r, ±fio) are approximately 8.3 and 2.9 (in units of Np)- 
The shaded profile shows the projection of N{r, ±fio) on the 
vertical plane. 

triguing to note that the impurity-induced STM signal 
is most strongly enhanced in the horizontal and vertical 
directions (the extrema directions of the rf-wave energy 
gap) at distances r Co from the impurity when the bias 
voltage passes through the resonance condition uj = ±Slo- 
The STM signal, proportional to N(r,ui) and shown in 
Fig. 2, has an unusually large maximum at r ~ Ip along 
these directions. Second, further away from the impurity 
~ ^o), the situation is reversed: the STM signal is 
enhanced along the diagonals for uj = — Slo (the hole-like 
resonance) and in the directions close to the diagonals 
for w = Slo (the electron-like resonance). For r ^ ^0, 
this enhancement should be large enough to be clearly 
observable, as illustrated in Fig. 3. However, at large 
distances, the impurity-induced contribution 6N(r, ±Slo) 
decays as making its detection experimentally de- 
manding. Third, particle-hole symmetry of the local 
density of states is manifestly broken by the impurity. 
Finally, our results are naturally understood by the ob- 
servation that, for c ^ 1, the problem effectively reduces 
into an isolated impurity and the superconductor with a 
vacancy. 



3 



0.0 1.0 3.0 3.0 4.0 



S; 0.2 







(a) ■ 




■=^- -1 - . 1 


1 






/In ' 1 

/ \ 
/ ^ 


(b) _ 




J \ \ 1 L. 





0.0 -1.0 -2.0 -3.0 -4.0 

Energy u/Q, 

FIG. 3. The local density of states N{r,Lo) of a rf-wave su- 
perconductor for (a) a; > and (b) a; < as a function of en- 
ergy Lo at the distance r/{o = 0.9 along the horizontal/vertical 
(solid lines) and diagonal (dashed lines) directions from the 
repulsive impurity (c = 0.1). The impurity resonance energy 
is Qo = 0.03Ao and ^o/ip = 6. 

This description leads to two impurity-induced fea- 
tures that simultaneously exist in the system: a state as- 
sociated with the impurity site itself and located at a very 
high energy (~ U), and a scattering resonance induced 
by the vacancy and describing a low-energy (SIq ^ Aq) 
"surface" state whose maximum weight is in the adja- 
cent sites to the defect. For non-zero c, these two fea- 
tures are weakly coupled rationalizing the resonance in 
N(r = 0,w) at w = — Slo- The resonant response of 
N(r 7^ 0,w) at w = ±Slo is the signature of the "surface" 
state, the virtual impurity-bound quasiparticle state. 

(ii) An anisotropic s-wave superconductor. As a pro- 
totypical example, consider the gap function A(k) = 
Ao|cos2i^|. Although it has the same bulk density of 
states as the rf-wave superconductor, there are no res- 
onance states regardless of the value of c [12]. This 
is a simplified version of the Anderson's theorem [1], 
which essentially states that disorder averages the energy 
gap over the Fermi surface; here the impurity scatter- 
ing leads to the same effect, resulting locally in an effec- 
tive isotropic energy gap of order Aq. The anisotropic 
s-wave superconductor behaves then like an ordinary 
isotropic s-wave superconductor regarding the resonance 
phenomenon. 

(Hi) A mixed s- and d-wave superconductor. Finally, 
consider a pairing state described by a gap function 
A(k) = Aq cos -|- As, where a small s-wave compo- 
nent is added to the rf-wave state, A^ ^ Aq. This 
type of admixture of two angular-momentum states is 
relevant in a rf-wave superconductor where the tetrago- 
nal symmetry is broken by an orthorhombic distortion. 
For \ui\ Aq, Gi(w) — As/Ao, and in the T-matrix 
equations c is replaced by an effective impurity strength 
2 = \J -\- (As/Ao)^. This implies that the strength of 
the impurity potential is effectively cut off at large values 



limiting its ability to generate a quasiparticle resonance. 
Note that Cj^\x,Lij) still approximately vanishes along 
the nodes of the energy gap, making the distorted rf-wave 
case similar to that of the tetragonal rf-superconductor 
so long as Aj/Aq < c. Thus, a small s-wave component 
would be invisible in an STM experiment. Controlling 
the degree of the orthorhombicity, the position and the 
sharpness of the resonance can be tunned. 

In summary, a central feature of low-energy quasiparti- 
cle states perturbed by an impurity is that they embody 
information about the underlying superconducting con- 
densate and its symmetry. An important consequence of 
a strongly-scattering impurity is the quasiparticle reso- 
nance which is sensitive to the phase of the pairing state 
and could be detected by STM. This property offers an 
intriguing possibility for the identification of the struc- 
ture of the pairing state in high-T^ superconductors. 
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